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q ' Abstract 

Effects of the mismatch between the hydrophobic length, d, of trans- 
membrane alpha helices of integral proteins and the hydrophobic thickness, 
■ Dh, of the membranes they span are studied theoretically utilizing a mi- 

croscopic model of lipids. In particular, we examine the dependence of the 
00 | period of a lamellar phase on the hydrophobic length and volume fraction 

of a rigid, integral, peptide. We find that the period decreases when a 
short peptide, such that d < Dh, is inserted. More surprising, we find that 
the period increases when a long peptide, such that d > D^, is inserted. 
The effect is due to the replacement of extensible lipid tails by rigid pep- 
tide. As the peptide length is increased, the lamellar period continues to 
increase, but at a slower rate, and can eventually decrease. The amount of 
£h [ peptide which fails to incorporate and span the membrane increases with 

the magnitude of the hydrophobic mismatch \d — Dh\. We explicate these 
behaviors which are all in accord with experiment. Predictions are made 
for the dependence of the tilt of a single trans-membrane alpha helix on 
^ ' hydrophobic mismatch and helix density. 

1 Introduction 



The interaction of integral proteins with the lipids in which they are embedded 
is of great importance for membrane function (Dumas et a!., 1999). One princi- 
pal governing this interaction is that the length of the hydrophobic segment of 
protein should closely match that of the membrane which it spans (Bloom and 
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Mouritsen, 1984). Among the evidence that this hydrophobic matching is used 
in membrane organization is that the various membranes between the Golgi and 
the plasma membrane have different thicknesses. Proteins can be routed through 
the secretory pathway by increasing their hydrophobic thickness via mutagene- 
sis and passing from one membrane to another more closely matching their new 
length (Pelham and Munro 93). The difference between the hydrophobic length 
of protein and membrane, denoted hydrophobic mismatch, affects inter alia, lat- 
eral segregation of proteins in membranes (Marsh 1995; Lehtonen and Kinnunen 
1997), the lipid melting transition (Piknova et al., 1993), and protein activity 
(Johannsson et al. 1981: Froud et al. 1986). Hydrophobic mismatch also affects 
the way in which the stabilty and the inclination of transmembrane helices change 
as functions of their hydrophobic length. Such information is very important in 
predicting transmembrane domains from potential protein sequences (Rost et al. 
1995; Edelman 1993), a topic becoming increasely important in Biology with the 
completion of the Human Genome Project. 

In an effort to elucidate such effects on a molecular level, Killian and co- 
workers (de Planque et al. 1998) investigated the effect of a series of hydrophobic 
peptides on the mean thickness of phosphatidylcholine membranes with differ- 
ent tail lengths. The peptides consisted of a sequence, whose length could be 
adjusted, of alternating leucine and alanine flanked on both sides by two trypto- 
phans. The latter prefer to reside just below the lipid head groups, and therefore 
serve as anchors for the peptide (Killian et al., 1996). The N- and C- terminii 
were blocked, e.g FmAW 2 (LA) n W 2 Etn. The results of this study which are of 
the most interest to us are as follows: 

• peptides whose hydrophobic thickness, d, is smaller than that of the hy- 
drophobic thickness, Dh, of the bilayer cause a reduction in the bilayer 
thickness. 

• peptides whose hydrophobic thickness is larger than that of the bilayer cause 
an increase in the bilayer thickness. As the peptides are made longer, the 
membranes continue to thicken, but the increase in the membrane thickness 
is always less than the peptide increment. 

• an increase in the magnitude of mismatch, \d — Dh\, whether due to peptides 
being too long or too short for the membrane, results in an increase in the 
fraction of peptides which fail to incorporate into the membrane. This effect 
has also been observed by Ren et al. (Ren et al. 1997, 1999). 

The first observation is easily understood on simple energetic arguments, but 
the second is not. One might have expected the longer peptide simply to insert 
at an angle such that the membrane thickness is undisturbed (Killian, 1998). 
However the thickness is disturbed, and grows with increased peptide length, but 
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does not track that increase identically. To understand this puzzling behavior 
and to isolate the various factors which bring about the net result of hydrophobic 
mismatch is the primary purpose of this paper. 

Many theoretical approaches have been applied to the effects of mismatch. 
Almost all of them are phenomenological, based on treating the membrane as an 
elastic sheet. Some have, however, included in these phenomenological descrip- 
tions some properties of the lipid tails, such as the ability of the tails to tilt to 
accommodate the perturbation of inserted proteins (May and Ben-Shaul, 1999). 
These approaches, which have recently been reviewed (May 2000), are valuable 
in clarifying several aspects of the general problem of membrane impurities, and 
have even been applied (May and Ben-Shaul, 1999) to the lamellar- inverted 
hexagonal transition induced by a sufficient concentration of peptide (van der 
Wei et al., 2000). Nevertheless, they lack a direct link to the molecular details of 
the system. A major theoretical advance was the work of Fattal and Ben-Shaul 
(Fattal and Ben-Shaul, 1993) who provided a molecular theory for the behav- 
ior of the lipid chains of the membrane. The peptides, with their hydrophobic 
length, were treated as providing a boundary condition on the configuration of 
the lipid chains. This molecular modelling was combined with phenomenological 
free energy contributions describing lipid headgroup repulsion and membrane sol- 
vent surface tension. In this paper we eschew phenomenological description and 
present a molecular theory which, from straight forward statistical mechanics, 
yields the free energy of the entire system, lipid and peptide. We utilize a molec- 
ular lipid model employed earlier (Li and Schick, 2000) and treat the peptide, 
which traverses the membrane as an alpha helix, as a rigid rod. We consider a 
lamellar phase formed by the lipid, and investigate the effect on the period of this 
phase due to the addition of peptide in small concentrations. As a consequence 
of our calculations, we are able to reproduce all of the above results, elucidate the 
reasons for the increase in membrane thickness when penetrated by long peptides, 
and delineate several different effect of hydrophobic mismatch. 

2 The Model and its Self Consistent Field So- 
lution 

We consider an anhydrous system of volume V consisting of lipids and peptides 
whose numbers, N\ and N p , are controlled by the fugacities z\ and z p , respectively. 
The lipids consist of a headgroup of volume t>h, and two equal- length tails each 
consisting of N segments of volume v t . Each lipid tail is characterized by a radius 
of gyration R g = (A^a 2 /6) 1 / 2 , with a the statistical segment length. The peptide 
consists of a rigid, hydrophobic, core of L segments each of length b and volume v c 
and a terminal group at each end of volume v e . We choose these end segments to 
be hydrophilic so that the peptide indeed models a hydrophobic segment within 
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an otherwise hydrophilic region. The hydrophobic length of the peptide is d = Lb. 

There are four local densities which specify the state of the system. We 
measure them all with respect to the convenient density v^ 1 . They are the 
number density of the lipid headgroup, 1( 3>h(r), and of the lipid tail segments, 
■y^ 1( 3?t(r), the number density of the peptide core segments, t>^ 1( l ) c (r), and the 
peptide end groups, t>^ 1( 3> e (r). Note that all number densities, $(r), are defined 
to be dimensionless. 

We consider repulsive contact interactions between these four elements. In 
the simplest case, the strengths of the repulsive interactions between the two hy- 
drophilic and the two hydrophobic elements are the same, kTv^x-, with k Boltz- 
mann's constant, and T the absolute temperature. We take the total energy of 
interaction to be 

/dv 
-[$ h (r) + 7e $ c (r)][ 7t $ t (r) + 7c $ c (r)]; (1) 

where we have introduced the relative volume of the tails, 7t = 2Nv t /vh, of the 
peptide cores, 7c = Lv c /vh, and of the peptide end groups 7c = v e /vh- Note 
that the interaction energies depend upon the local volume fractions 7e $ e etc. as 
opposed to the local number densities $ e etc. (Williams and Fredrickson, 1992). 

In addition to this local repulsion, we include the hard core interactions be- 
tween all particles in an approximate way by imposition of a local incompressibil- 
ity constraint, i.e., that the sum of the volume fractions of all components must 
be unity everywhere: 

A(r) = $ h (r) + 7t $ t (r) + 7c $ c (r) + 7c $ c (r) - 1 

= 0. (2) 

As shown earlier (Li and Schick, 2000), the partition function of the system can be 
written in the form in which the four fluctuating densities, instead of interacting 
directly with one another, interact indirectly via four fluctuating fields, here 
denoted W h , W t , W e , and W c . Self consistent field theory results when the 
fluctuating fields and densities are approximated by those values which extremize 
the free energy, Q, of the system in the presence of these fields. This free energy 
has the form 

n Qi[W h ,W t ] Q P [W C ,W C ] v h 

Si — — Z\ z n \- tli 

kTV V p V kTV 

/dv 
-[W h (r)$ h (r) + W t (r)$ t (r) + W c (r)<S> c (r) + W c (r)$ c (r)] 

- /*S(r)A(r). (3) 

Here Qi[Wh, W t ] is the partition function of a single lipid in external fields Wh, and 
W t , and Q p [We, W c ] is the partition function of a single peptide in the external 
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fields W e and W c . Note that a Lagrange multiplier S(r) has been introduced 
to enforce the incompressibility constraint of Eq. |2|. The functions W&, $h etc. 
which extremize this free energy will be denoted by their corresponding lower 
case letters and 0h- They are obtained from the five equations for the fields 

= X[7t0t(r)+7c0c(r)]+£(r), (4) 

w t (r)/7t = x[0 h (r) +7 c c (r)] + £(r), (5) 

^ c (r)/7e = x[7t0t(r) +7c0c(r)] +^(r), (6) 

^c(r)/7c = x[0 h (r)+7c0c(r)]+e(r), (7) 

1 = h (r)+7 t t (r) + 7c c (r)+7 c c (r). (8) 

The field £ can be easily eliminated, while Eqs. © and ([?]) imply w t (r)/j t = 
u> c (r)/7 c , so that one deals essentially with three equations. The four densities 
are all functionals of the above fields except £ and, therefore, close the cycle of 
self-consistent equations: 

, ( v n SQi[w h ,w t ] 

(p h {r)[w h ,Wt\ = -Zi—z — — , (9) 

<k(r)[w h ,w t ] = -z x , 10 

5iut(r) 

c (r)[w e ,w c J = -z p ^ w 73 — > I 11 ) 

C r K ? w c = ~*p — 7 — 7-v — • 12 

ow c (r) 

With the aid of the above equations, the self consistent, or mean field, free energy 
f2 m f, which is the free energy function of Eq. [5] evaluated at the self consistent 
field values of the densities and fields, can be put in the form 

kT 

-il mi (z h z p ,T) = — (ziQi[w h ,w t ] + z p Q p [w e ,w c ] + V£ )+E[(j>h, fa, <f> e , C ] (13) 

where we have defined V£o = / £( r )<^ r - All the fields, w, £ > and densities, 0, are 
functions of the activities, z\, z p and temperature. All of the above is a simple 
extension of previous procedure (Li and Schick, 2000). 

There remains only to specify how the partition functions of the lipids and of 
the peptides are calculated. One defines the end-segment distribution function 
of the lipid q(r, s). Because the lipid tails are treated as completely flexible, this 
function satisfies the modified diffusion equation 

- 2^v 2 g(r, s) + [w h (v)5(s - 1/2) + w t (v)} q(v, s) = 0, (14) 
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with initial condition 

g(r,0) = l. (15) 
From this function, one obtains the partition functions of the lipids, 

Qx = J dr g(r,l), (16) 

and, from Eqs. @ and (|10D, the head and tail densities 

h (r) = zi exp{-w h (r)}g (r, q (r, --^ , (17) 

t (r) = zi [ ds q(r, s)q(r, 1 - s) (18) 
Jo 

To obtain the partition function of the peptide, one defines its end-segment distri- 
bution function q p (r, n, s), where n is a unit vector which specifies the orientation 
of the peptide. Because the peptide is rigid, and of length Lb, this function sat- 
isfies the equation (Wang and Warner, 1986) 

9gp( ^ n,g) + Lbh • V r g p (r, n, s) + {w c (v)[6(s) + S(s - 1)] + ^ c (r)}g p (r, n, s) = 0, 

(19) 

with initial condition 

g p (r,n,0 + ) = exp[-w (r)]. (20) 
From this function, one obtains the partition function of the peptide, 

Q p = JdrJ dng p (r,n,l), (21) 

and, by means of Eqs. (|TTD and ([T2]), its end and core densities, 

<f) c (r) = 2z p Jdhq v (r,h,l), (22) 

,1/2 , 

c (r) = 2z p ds J dh q p (r, h, s)q p (r, -h, 1 - s). (23) 

To summarize: there are five self-consistent equations to be solved for the five 
fields Wh(r), ti>t( r ), w c (r), w c (r), and £(r). They are Eqs. (|[) to (H). The fields 
depend on the four densities 0h( r ), </>t( r ), 0e( r ), and <t> c { r ), which depend, in turn, 
on these fields. The densities are given by Eqs. [T7|, [18], and |2^. Once the 
fields and densities are obtained, the free energy follows from Eq. |l^. 

We are interested in the way in which the peptides affect, on the average, 
a periodic array of lipid bilayers, that is, a lipid lamellar phase. We therefore 
assume that the lamellae are uniform in their plane, and vary only in the normal 
direction, specified by the coordinate z. 
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In the limit of vanishing peptide density, it is sufficient to solve for the fields 
and densities of the pure lipid bilayer, and then to solve Eq. ( |19"|) for the peptide 
end-segment distribution function in the presence of those densities. This is easily 
done in real space, since the solution of ([TJ) is just 

g p (r, n, s) = exp — / dt w c (r + th)[5(t) + 5(t — 1)] + w c (r + th) , 
L Jo 

and the fields are provided by Eqs. @ and (0) with <p e = (f> c = 0. 

At non-zero peptide densities, the full set of self-consistent equations must be 
solved, and it is more convenient to expand all functions of the position r and 
the director n in terms of a complete set of basis functions (Wang and Warner 
1986, Matsen 1996). 

oo 

9( T ,^) = J2 9i,mfi,m(Zi cos 0), (24) 

l,m=0 

where cos 9 is the projection of the unit vector n onto the z axis. A convenient 
set is 



/,,„(*, cos 0) = (2J + l) 1/2 fl(cos0), 

'27imz\ 



f hm {z,cos9) = (4l + 2) 1 / 2 cos(- 



P/(cos I 



f ljm (z,cos9) = {41 + 2) 1/2 sin (-^p) ^ (cos 9) odd m, (25) 

where Pi is the /'th Legendre polynomial, and D is the period of the lamellae. 
The latter is determined by minimization of the free energy with respect to it. 
Details of the procedure for solving the self-consistent equations in this basis can 
be found elsewhere (Matsen 1996). Of importance here is that the expansion into 
the infinite set of basis functions must be truncated to an expansion in a finite 
number of such functions. We have utilized 30 values of / and 10 values of m, 
or a total of 300 functions. The errors in the free energy brought about by this 
truncation are less than 0.1%. 

The parameters we have chosen for our calculations are as follows. The lipid 
is characterized by the volume of the headgroup, = 370A 3 , and the volume of 
the tails, 2Nv t = 999A 3 . For comparison, the volume of the phosphatidylcholine 
headgroup is 337A 3 and that of two tails with seventeen carbons and one double 
bond each is 985A 3 (Armen at al., 1998). The radius of gyration of the tails was 
taken to be R g = 12. 3A which was found to be reasonable in a previous study 
(Li and Schick 2000). The peptide is characterized by the volume of the end 
groups, which we took to be f e = 555A 3 each, the volume of each of its core 
units, v c = 96. 2A 3 , and the length of each core group b which we took to be the 
length of each amino acid in an alpha helix conformation, 1.5 A. For comparison, 
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the volume of the two tryptophans and two alanines, which were only a portion 
of the end groups used by de Planque et al. (de Planque et al., 1998), is 460A 3 , 
while the average of the volumes of the core alanine and leucine units is 95. 5A 3 . 
The interaction strength between hydrophilic and hydrophobic elements is such 
that X7t = 20. We have taken the number of peptide units to vary from L = 10 
to 40. 

3 Results 

We first consider the limit in which the density of peptide is vanishingly small. 
The volume fraction distribution of the lipid headgroups and tails in the anhy- 
drous lipid bilayer are shown in Fig. 1. The period D, in units of the radius of 
gyration of the lipid tails, R g is D/R g = 2.831. The thickness of the hydrophobic 
region, D^, as defined by those points at which the volume fraction of the tails 
7t$ t = 0.5, is D h /D = 0.757. 

It is of interest in this limit to determine whether the inserted peptide spans 
the bilayer and, if so, whether it inserts normal to the bilayer or at an angle, 9, to 
it. This is readily determined. We compute the probability distribution (Matsen, 
1996) 

Pl(z, cos 9) = q p (z, cos 9, 1)/ J dzd(cos9)q p (z,cos9,l) 

of peptide ends which are at an angle 9 with the bilayer normal, when that end 
is located at position z within the bilayer. The coordinates z/D = and 1 cor- 
respond to the center of the tail region, as in Fig. 1. The probability distribution 
can be calculated for peptides of different length Lb. It is shown in Fig. 2 for 
an L = 20 peptide, for which Lb/D h = 1.143, that is, somewhat longer than the 
hydrophobic thickness of the bilayer. One sees two major peaks, both of which 
correspond to a peptide which spans the bilayer. One corresponds to a peptide 
inserted almost normally, at an angle such that cos 9 pa —0.9 and with one end at 
z/D pa 0.4 This peptide would pass through the tail region at z/D = 0. Ends of 
peptides in the adjacent lamellae near z/D ~ 0.6 are characterized by an angle 
7r — 9 so that the value of the cosine is the negative of that of the first peak. Thus 
the second peak simply describes the other end of the rigid, (and periodically 
repeated), peptide. 

In order to illuminate the behavior of this angle of insertion, we integrate the 
probability distribution of Fig. 2 over all spatial positions for a fixed angle to 
obtain the probability distribution P^(cos^). It is shown in Fig. 3 for several 
values of the hydrophobic length of the peptide Lb divided by the hydrophobic 
thickness of the bilayer, D^. From this distribution we obtain the average angle of 
peptide insertion and the most probable angle of insertion as a function of Lb/D h . 
These quantities are plotted in Fig. 4 in dashed and solid lines respectively. We 
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observe that peptides with hydrophobic lengths (Lb/D^) < (L*b/D\/) = 1.07, 
insert normal to the bilayer. One might have expected this ratio to be unity, 
but it must be recalled that our definition of the hydrophobic thickness of the 
membrane, Dh, in terms of equal head and tail volume fractions is a somewhat 
arbitrary one. Peptides which are longer than L*b insert an a non-zero angle to 
the membrane normal. From energetic arguments alone, one might expect that 
cos 9 oc 1/L for L > L*. Indeed in our calculation, one sees the dependence on 
Lcos6 in the second term of Eq. (Jl9|). However the incompressibility condition 
depends only on the volume of the peptide, proportional to L, not on its orien- 
tation. Thus the behavior cos 9 oc 1/L should only be a simple approximation 
to the actual behavior. That this is indeed so is seen in Fig. 4 where we have 
plotted this dependence as a dotted line. One sees that our results deviate from 
this simple description. 

Returning to Fig. 2, we also observe two smaller peaks at cos 9 = corre- 
sponding to a fraction of peptides which do not traverse the lamellae, but are 
parallel to it. As our model peptide is predominantly, but not completely, hy- 
drophobic, the non-traversing peptides are found somewhat below the head-tail 
interface in the tail region of the bilayer. In order to determine how this amount of 
non-traversing peptides varies with the peptide length, we calculate the fraction 
of peptides with an end at an angle < cos6> < 0.2, 



This fraction is plotted vs. Lb/D^ in the inset of Fig. 3. We see that the 
fraction which do not insert across the membrane increases with the mismatch 
between peptide and bilayer hydrophobic thicknesses, (Lb/Dh) — 1, irrespective 
of the sign of the mismatch. The largest fraction of inserted peptides occurs for 



We now consider non-zero peptide densities. We have calculated the period, 
D, of the lamellar phase as a function of peptide volume fraction, x p , for values 
of x p < 0.04 and find that D varies essentially linearly with it; i.e. 



where -D(O) is the period of the lamellar phase in the absence of peptide. In Fig. 
5 we plot the rate of bilayer thickening, R, versus the relative peptide thickness 
Lb/Dh- We see that the insertion of short peptides, Lb/Dh < 1.17, causes the 
period, and therefore the bilayer thickness, to contract, while the insertion of 
long peptides causes it to increase. As the peptides become longer, the rate of 
increase becomes smaller. Eventually the insertion of sufficiently long peptides, 
Lb/Dh > 1-6, causes the bilayer thickness to decrease. 




(26) 



Lb/D h « 1.15. 



D(x p )/D(Q) ~ 1+Rx, 



(27) 
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4 Discussion 



Our results replicate and illuminate the experimental results listed in the Intro- 
duction. We recall and discuss them in turn. 

The insertion of peptides whose hydrophobic thickness is less than 
that of the bilayer causes a reduction in the bilayer thickness. 

Our calculation reproduces this result. The behavior is clearly due to simple 
energetic considerations. 

The insertion of peptides whose hydrophobic thickness is greater 
than that of the bilayer causes an increase in the bilayer thickness. As 
the peptides are made longer by a certain amount, 5d, the membranes 
continue to thicken, but the increase in membrane thickness is less 
than 8d. 

We reproduce these results. As noted earlier, their origin is not obvious, for 
energetically the peptides would be satisfied to insert at an angle such that the 
bilayer would not deviate from its original thickness at all. The reason for this 
behavior, therefore, is clearly entropic. To understand it, we need only recall the 
reason the bilayer thickness takes the value it does in the absence of peptide. 
The repulsive interaction between lipid headgroups and tails tends to crowd the 
head groups together which, from the constraint of incompressibility, causes the 
lipid tails to stretch. This tendency is opposed by the loss of tail entropy such 
stretching brings about. The equilibrium thickness of the bilayer results from 
a balance of these two tendencies. Thus in the equilibrium configuration, the 
lipid tails are stretched (Fattal and Ben-Shaul, 1993). Peptides, being rigid, 
do not stretch, and therefore lose no such entropy if the bilayer thickens, while 
they displace lipid tails which restrain the membrane from thickening. Thus it 
can be understood that the insertion of longer peptides causes the width of the 
bilayer to increase. Within this mechanism, however, there is no reason that 
a given increase in the length of the peptide should result in a corresponding 
increase in the thickness of the bilayer. The effect of replacing elastic lipids with 
rigid peptides predicts that the insertion of a much bulkier rigid object, such as 
gramacidin, will have a larger effect in increasing the bilayer thickness than will a 
less bulky one, such as a simple peptide. This is again in accord with experiment 
(de Planque, 1999). 

An increase in the magnitude of mismatch, whether due to peptides 
being too long or too short for the membrane, results in an increase 
in the fraction of peptides which fail to incorporate in the membrane. 

Our calculation reproduces this as shown in the inset to Fig. 3. Short peptides 
tend not to insert for energetic reasons. Long peptides tend not to insert even 
though there is no energetic penalty to do so. Presumably, they fail to insert due 
to the loss of lipid tail entropy, which increases with the length of the hydrophobic 
portion of the inserted peptide. 
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In addition to reproducing these experimental results which demonstrates its 
efficacy, our model also yields predictions. As noted earlier, we found that inser- 
tion of peptides whose hydrophobic length is greater than the hydrophobic thick- 
ness of the membrane causes the membrane thickness to increase. As the peptide 
is made longer, the rate of increase of membrane thickness with peptide length, 
denoted R in Eq. (0), decreases in accord with experiment (de Planque,2000). 
Strikingly, our model predicts that insertion of peptides which are very long com- 
pared to the hydrophobic thickness of the membrane (Lb/Dh > 1.6 in Fig. 5) 
actually causes the membrane thickness to decrease. We believe the reason for 
this is that a significant fraction of these long peptides do not traverse the mem- 
brane, but lie parallel to it. Because the peptides have a much larger end group 
relative to its core than do the lipids, they create relatively more free volume for 
the lipid tails to fill, i.e. they effectively increase the area per lipid head group. 
As the system is incompressible, this effect tends to make the bilayer thickness 
decrease. We have recently learned that this unusal thinning of the bilayer on 
the addition of relatively long peptides has been observed (Killian, 2001). 

Finally, our model makes predictions about the insertion of peptides at very 
low densities. Short peptides, whose hydrophobic thickness is less than the hy- 
drophobic thickness of the bilayer, insert normal to it. Peptides with a hydropho- 
bic thickness greater than that of the bilayer insert at an angle to the normal 
which grows with the hydrophobic mismatch. We expect the same behavior for 
proteins which span the membrane with a single alpha helix. As we have ignored 
local effects in the plane of the membrane, the actual angle of insertion will dif- 
fer somewhat from that which we have calculated, but the qualitative behavior 
will not be changed. Except for very large mismatches, our theory predicts the 
tilt to decrease with increasing peptide concentration. Although there is much 
data on helix tilt in specific systems, there appears to be no attempt at a sys- 
tematic correlation of it with hydrophobic mismatch. Such data would be most 
interesting. 
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thesis of M. R. R. de Planque. This work is supported by the NSF under grant 
DMR9876864. D.D. acknowledges support from a Spanish Ministry of Education, 
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11 



REFERENCES 

Armen, R., O. Uitto, and S. Feller. 1998. Phospholipid component volumes: 
determination and application to bilayer structure calculations. Biophys. 
J. 75:734-744. 

de Planque,M. R. R. 2000. Hydrophobic matching and interfacial anchoring of 
transmembrane peptides as determinants for membrane organization. Ph. 
D. Thesis, University of Utrecht. 

de Planque, M. R. R., D. V. Greathouse, R. E. Koeppe II, H. Schaffer, D. Marsh, 
and J. A. Killian, 1998. Influence of lipid/peptide hydrophobic mismatch on 
the thickness of diacylphosphatidylcholine bilayers. A 2 H NMR and ESR 
study using designed transmembrane a-helical peptides and gramacidin A. 
Biochemistry 37:9333-9345. 

de Planque, M. R. R. , J. A. W. Kruijtzer, R. M. J. Liskamp, D. Maarsh, D. 
V. Greathouse, R. E. Koeppe II, B. de Kruijff, and J. A. Killian, Differ- 
ent membrane anchoring positions of tryptophan and lysine in synthetic 
transmembrane a-helical peptides. J. Biol. Chem. 274:20839-20846. 

Dumas, F., M. C. Lebrun, and J.-F. Tocanne. 1999 Is the protein/lipid hy- 
drophobic matching principle relevant to membrane organization and func- 
tions? FEBS Letters 458:271-277. 

Edelman, J. 1993. Quadratic Minimization of Predictors for Protein Secondary 
Structure: Application to Transmembrane a- Helices. J. Mol. Biol. 232:165- 
191. 

Fattal, D. R. and A. Ben-Shaul. 1993. A molecular model for lipid-protein 
interaction in membranes: the role of hydrophobic mismatch. Biophys. J. 
65:1795-1809. 

Froud, R.J., C. R. A. Earl, J. M. East, and A. G. Lee. 1986. Effects of lipid fatty 
acyl chain structure on the activity of the (Ca2+- Mg2+)-ATPase. Biochim. 
Biophys. Acta 860:416-421. 

Johannsson, A., G. A. Smith, and J. C. Metcalfe. 1981. The effect of bilayer 
thickness on the activity of (Na+ + K+)-ATPase. Biochim. Biophys. Acta 
641:416-421. 

Killian, J. A., I. Salemink, M.R.R. de Planque, G. Lindblom, R.E. Koeppe 
II., and D.V. Greathouse. 1996 Induction of nonbilayer structures in dia- 
cylphosphatidylcholine model membranes by transmembrane a-helical pep- 
tides: Importance of hydrophobic mismatch and proposed role of trypto- 
phans. Biochemistry 35:1037-1045. 



12 



Killian, J. A. 1998 Hydrophobic mismatch between proteins and lipids in mem- 
branes. Biochim. Biophys. Acta 1376:401-416. 

Killian, J. A. 2001., private communication. 

Lehtonen, J. Y. A. and P. K. J. Kinnunen. 1997. Evidence for phospholipid 
microdomain formation in liquid crystalline liposomes reconstituted with 
Esherichia coil lactose permease. Biophys. J. 72:1247-1257. 

Li, X.-J., and M. Schick. 2000. Theory of Lipid Polymorphism: Application to 
Phosphatidylethanolamine and Phosphatidylserine. Biophys. J 78:34-46. 

Marsh, D. 1995. Lipid-protein interactions and heterogeneous lipid distribution 
in membranes. Mol. Membr. Biol. 12:59-64. 

Matsen, M. W. 1996. Melts of semiflexible diblock copolymer, J. Chem. Phys. 
104:7758-7764. 

May, S. 2000. Theories of structural perturbations of lipid bilayers. Curr. 
Opinion in Coll. and Interface Sci. 5:244-249. 

May, S. and A. Ben-Shaul. 1999. Molecular theory of lipid-protein interaction 
and the L a — Hn transition. Biophys. J. 76:751-767. 

Mouritsen, O. G., and M. Bloom. 1984. Matress model of lipid-protein interac- 
tions in membranes. Biophys. J. 46:141-153. 

Pelham, H.R. and S. Munro. 1993. Sorting of membrane proteins in the secre- 
tory pathway. Cell 75:603-605 

Piknova, B. E. Perochon, and J. F. Tocanne. 1993. Hydrophobic mismatch and 
long/range protein lipid interactions in bacteriorhodopsin/phosphatidylcholine 
vesicles. Eur. J. Biochem. 218:385-396. 

Ren, J., S. Lew, Z. Wang, and E. London. 1997. Transmembrane orientation of 
hydrophobic a-helices is regulated both by the relationship of helix loength 
to bilayer thickness and by the cholesterol concentration. Biochemistry 
36:10213-10220. 

Ren, J., S. Lew, J. Wang, and E. London. 1999. Control of the transmembrane 
orientation and interhelical interactions within membranes by hydrophobic 
helix length. Biochemistry 38:5905-5912. 

Rost, B., R. Casadio, P. Fariselli, and C. Sander. 1995. Transmembrane helices 
predicted at 95% accuracy. Protein Sci. 4:521-533. 



13 



van der Wei, P. C. A., T. Pott, S. Morein, D. V. Greathouse, R. E. Koeppe 
II, and J. A. Killian. 2000. Tryptophan-anchored transmembrane pep- 
tides promote formation of nonlamellar phases in phosphatidyletanolamine 
model membranes in a mismatch-dependent manner. Biochemistry 39:3124- 
3133. 

Wang, X. and M. Warner. 1986. Theory of nematic backbone polymer phases 
and conformations. J. Phys. A:Math. Gen. 19:2215-2227. 



14 



Figure 1 Volume fraction distribution of the lipid headgroups and tails in the 
lamellar phase plotted vs. coordinate z perpendicular to the lamellae. The 
wavelength of the phase is D. z/D = 0.5 corresponds to the center of the 
head region, and z/D = and 1, to the center of the tail region. The 
interaction strength between hydrophobic and hydrophilic entities is such 
that X7t = 20. 

Figure 2 Normalized probability distribution, P L (z,cos9) of the orientation, 
cos 9, and location, z, of the peptide ends for a L = 20 peptide, of length 
Lb/Dh = 1.143, with D^ the hydrophobic thickness of the lamellae. The 
center of the lipid tail region is at z/D = and 1. Contour values are given 
in the legend. 

Figure 3 Normalized probability distribution Pj,(cos#) of the orientation of the 
peptide ends vs. cos 9 for various relative peptide lengths Lb/D h , as noted 
on each curve. The lengths correspond to L = 16 to L = 26 in increments 
of 2, plus two extreme cases, L — 10 and L = 40. The inset shows the 
fraction, J, of peptides which do not insert into the bilayer as a function of 
Lb/D h . 

Figure 4 Full line: location of the peak in the probability distribution function 
P(cos8) which occurs at non-zero values of cos#; dashed line: mean value 
of cos#; dotted line: simple 1/L behavior. Curves are given as functions of 
peptide length in units of D^ (lower x scale) and L (upper x scale). 

Figure 5 Rate, R, of increase of relative lamellar period with volume fraction, 
x p , of peptide vs. relative hydrophobic length of peptide. 
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